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We examine the taste structure of eigenvectors of the staggered-fermion Dirac operator. We 
derive a set of conditions on the eigenvectors of modes with small eigenvalues (near-zero modes), 
such that staggered fermions reproduce the 't Hooft vertex in the continuum limit. We also show 
that, assuming these conditions, the correlators of flavor-singlet mesons are free of contributions 
singular in 1/m, where m is the quark mass. This conclusion holds also when a single flavor of 
sea quark is represented by the fourth root of the staggered-fermion determinant. We then test 
numerically, using the highly improved staggered-quark action, whether these conditions hold on 
realistic lattice gauge fields. We find that the needed structure does indeed emerge. 
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I. INTRODUCTION 

Lattice QCD has made several notable strides over 
the past decade. A wide variety of calculations with 
2-\-l flavors of sea quarks (corresponding to up, down, 
and strange) have been found to agree with experimen- 
tal measurements within 2% [1]. Charmed- meson de- 
cay constants j2J, semileptonic form factors [3], and the 
masses of the Be [1] and rji, [S] mesons have been com- 
puted before being confirmed by measurements from ex- 
periments. Calculations at nonzero temperature have 
shown that QCD possesses not a first-order phase tran- 
sition but a smooth crossover 6^, with implications for 
heavy-ion collisions and a cooling universe. Some of the 
most precise determinations of the strong coupling [7] , 
quark masses |5], and flavor-changing couplings [HI ITU] 
come from lattice QCD. It is impractical to cite every 
development here, but recent reviews [HI [12] cover the 
breadth of progress well. 

The results listed above pifTU] have been obtained us- 
ing staggered fermions [T31 [TJ] for the sea quarks, be- 
cause this approach is numerically the fastest [15j . In the 
continuum limit, one staggered-fermion field yields four 
species with a quantum number nowadays called "taste." 
In numerical lattice gauge theory, sea quarks are repre- 
sented by a determinant, for staggered fermions, 



Det (D. 

4 



stag 



(1) 



where Dstag denotes the lattice Dirac operator (see be- 
low), TO is the quark mass, and the subscript 4 is a re- 
minder that the natural outcome is 4 tastes. To simulate 



*Email: g.donald@physics.gla.a c.uk | 

tEmail: c.davies@physics.gla.ac.uk 

tEmail: efollana@unizar.es 

§ Email: as k@fnal.goV| 



a single species of given mass with staggered fermions, 
the (4-taste) determinant representing the sea is replaced 

with [in] 



Det [D; 

. 4 



stag 
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1/4 



(2) 



Below we shall refer to the systems using ([T]) and ^ 
as "unrooted" and "rooted" staggered fermions, respec- 
tively. As far as we know, there is no controversy that 
lattice gauge theory with unrooted staggered fermions ([T]) 
defines a four-species continuum gauge theory. 

The fourth root is controversial, however, because it 
is not standard quantum field theory. The arguments 
supporting its validity hinge on structural properties of 
unrooted staggered fermions, which suggest that the con- 
tinuum limit of Det4(L'stag + m) in Eq. ^ factors into 
four equivalent determinants [T7U20j . This factorization 
is verified in weak-coupling perturbation theory, where 
the 1/4 from the exponent multiplies each fermion loop. 
Weak coupling also suggests how the symmetries of four 
species emerge in the continuum limit. In simplified but 
similar systems where one can retain analytical control, 
the rooted determinant is valid [211123) . Extensive nu- 
merical studies elaborate how the procedure works in the 
Schwinger model ^24j 25 . Straightforward analysis of the 
hadron mass spectrum as a function of lattice spacing 
and quark mass, using chiral perturbation theory, sub- 
stantiates this picture in detail [111 I26j . Further nonper- 
turbative evidence comes from studying the eigenvalues 
of the staggered-fermion operator Dgtag: demonstrating 
that they appear in nearly degenerate quartets [?7H50] . 
On lattice gauge fields with nonzero topological charge, 
sets of quartets with eigenvalues near zero emerge. The 
number of quartets and their chirality satisfy the index 
theorem f^l5D]. 

One issue that has not been fully addressed is the be- 
havior of flavor-singlet mesons. Direct calculations of the 
flavor-singlet meson masses are difficult [3TH51], because 



2 




tex in the continuum limit, we find that the eigenvec- 
tors must satisfy the same properties derived in Ref. [19] , 



isospin 1 



isospin 



FIG. 1: Pattern of flavor and taste quantum numbers in 
the (pseudoscalar) meson spectrum with two flavors and four 
tastes. The flavor-nonsinglet (isospin 1) mesons are split by 
small lattice artifacts. The flavor-singlet (isospin 0) taste- 
nonsmglet mesons are no different. The flavor-taste singlet, 
however, receives a contribution from mixing with purely glu- 
onic states, an effect studied in Refs. 1311 1321. 



they entail a contribution in which the quark-antiquark of 
the meson annihilates into giuons, and the gluons recre- 
ate the quark-antiquark pair. With staggered fermions, 
it is crucial to bear in mind that only the flavor-taste 
singlet can undergo this process. Low-energy gluons are 
taste singlets, so a meson with nontrivial taste cannot 
annihilate into them. The spectrum with two flavors is 
sketched in Fig. [TJ illustrating the roles of the flavor and 
taste quantum numbers. 

Building on the eigenvalue studies, this paper addres- 
ses a specific concern, namely that flavor-taste-singlet 
correlators could diverge as a power of to as m — 0. 
Such behavior would be a clear failure of rooted staggered 
fermions. We find fault with key steps in an attempted 
derivation of this claim [551 [55] , which uses the 't Hooft 
vertex [371 EE] to try to understand the role of near-zero 
modes. A complementary examination of the same cor- 
relators reduces the problem to certain properties of the 
near-zero modes' eigenvectors [19] . Then contributions 
from connected and disconnected correlators cancel the 
divergent behavior; with the correct combinatoric fac- 
tors [T51 dni [SHHII] ) the cancellation holds even with the 
rooted determinant of Eq. ([2| . 

In this paper, we derive the staggered-fermion 't Hooft 
vertex directly from the functional integral, both for un- 
rooted and rooted staggered fermions. If unrooted stag- 
gered fermions are to obtain a four-species 't Hooft ver- 



namely Eqs. (31 ) and (32 1 below. References [JDJUI] tac- 
itly assumed these properties, but we examine the eigen- 
vectors numerically, plotting the quantities that enter the 
't Hooft vertex and the fiavor-taste-singlet correlators. 
We find that they behave in precisely the way needed for 
unrooted and rooted staggered fermions to yield four or, 
respectively, one species in the continuum limit. 

The rest of this paper is organized as follows. Section[TT| 
discusses staggered fermions and some of the complaints 
and concerns about Eq. Section III reviews the con- 
tinuum 't Hooft vertex and its symmetries, constructs 
the staggered-fermion 't Hooft vertex, and sets up the 
problem of fiavor-taste singlets. This discussion also pin- 
points where the analysis of Refs. [351 goes astray. 



Section [TV] explains details of our numerical setup, gives 
our lattice results, and discusses their implications. The 
data speak for themselves: they clearly show that the 
needed structure emerges dynamically, ever more so for 
smaller lattice spacing. Section \V\ gives our conclusions. 
It seems to us that the rooted staggered sea has passed 
another test in its usual way of relying on properties of 
the unrooted theory. Appendix [Xj contains some cum- 
bersome notation that lends technical completeness to 
Sees. HH |III[ and |IV| Appendix [B] writes out improved 
actions explicitly. Appendix [Cj remarks on issues of sec- 
ondary importance, raised in Refs. [351155] . 



II. STAGGERED FERMIONS 

In this section, we review unrooted staggered fermions, 
because the way that four species emerge is central to 
any argument that Eq. ([2]) is a valid regulator for one 
species. We are careful to distinguish between flavor and 
taste; the former is a label decoupled from the gauge 
interaction; the latter is a property of staggered fermions, 
described below. 

Below we use improved actions to check numerically 
whether the dynamics of staggered fermions are as ex- 
pected. For the discussion here, it is enough to start 
with the original, unimproved lattice action (14j: 



ia^ X! ^t^i^)x{x) [U^,{x)x{x + jia) - Ul{x - ila)x{x - fia)] + ma^ ^ x{x)x{x), 



(3) 



r 



where a is the lattice spacing, xi^) ^^'^ x{^) ^-re gauge- 
group multiplets of Grassmann numbers for lattice site x, 
U^{x) is a lattice gauge held connecting sites x and x+fla 
(such that S'stag is gauge invariant), to is the bare mass, 
/t is a unit vector in the /i direction, and /i G {1, 2, 3, 4}. 



The staggered-fermion fields carry no Dirac index, and 
sign factors appear instead of Dirac matrices: 



p < it p / 



(4) 
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The staggered Dirac operator -Dgtag is defined by writing 
x,y 

The determinant (fTl) follows from integrating the func- 
tional integral over (X: x)- 

"S'stag is invariant under shifts: 



where 



Xi^) Cf^{x)x{x + fia) 
Ui/ (x) I— ?> U^{x + fia) \lv 



(6) 



(7) 



Acting on fermion fields, SySy, 



-SuSu- This built-in 



Clifford group r4 is the origin of the four species in the 
continuum limit and their quantum number taste. Acting 
on gauge fields, SyS^,^ — +S^Si^, from which it follows 
that low-momentum gauge fields are taste singlets. With 
Uf fiavors of (x,x) — so Anj species in all — there is still 
only one gauge field and, thus, only one 

The kinetic term (for n f flavors) is also invariant under 
a U(n/) symmetry group. 



X{x) 1-^- x{x)e' 



<p"T°e(a;) ! 



(8) 



where the are (anti-Hermitian) flavor generators, in- 
cluding flavor-singlet T" = 



and 



(9) 



Three crucial properties of these symmetries ([8]) are that 

1. they are exact even at nonzero lattice spacing a; 

2. they are nonsinglets with respect to taste; 

3. they imply that the eigenvalue spectrum of -Dstag 
is pure imaginary and symmetric about 0. 

The first property means that these symmetries cannot 
be anomalous, so they cannot be germane to the index 
theorem. The second property means that the lack of 
anomaly is good: in QCD, species-nonsinglet symmetries 
do not have anomalies. A corollary of the third property 
ensures that if iX is an eigenvalue of Dstag with eigen- 
vector f{x), then ~iX is also an eigenvalue, now with 
eigenvector e{x)f{x). This corollary plays an important 
role in Sec. |IV[ Unfortunately, the connection between 
property [3] and the spectrum sometimes leads, it seems, 
to the fiavor-singlet Ue being misidentified as the ana- 
log of continuum QCD's anomalous 11^(1)- The first two 
properties mean, however, that even the flavor-singlet 
cannot be related to U^(l). 

The analog of the U^(l) is a flavor and taste singlet. 
It is explicitly broken for a ^ but restored — apart 
from the anomaly and mass terms — as a — > |14j . This 



mechanism is familiar in lattice gauge theory the 
same happens with Wilson fermions |43]. As a — 0, an 
anomalous Ward identity emerges with axial-vector cur- 
rent, Aj{x), and pseudoscalar density, Pi{x), that are 
taste-flavor singlets [HI HI] . The subscript / denotes the 
trivial representation of the shift symmetries Q, also 
called the taste-singlet representation. 

The way flavor-taste symmetries emerge is crucial to 
the validity of staggered fermions. In particular, 

r4 X SUviuf) C SUy(4n/), (10) 
Ue(n/)^U(n/)®C'cSU^(4n/), (11) 

where the symmetries on the left are exact (or softly 
broken) for S'stag, and those on the right are desired for 
continuum QCD. The SUy(n/) on the left-hand side of 
Eq. ( 10 1 is the obvious flavor- number symmetry of Eq. ^ 
for Uf fiavors of equal mass. The on the left-hand side 
of Eq. (ITT]) denotes the taste- nonsinglet nature of U^. 



The pattern of symmetry appears most vividly, both 
for nonzero a and as a — ^ 0, in the meson-mass spec- 
trum. Meson operators can be written X-T^Xi where 
denotes various choices of sign factors rj and paral- 
lel transport within a hypercube, such that the bilinear 
transforms under the bosonic representation of the 
shift symmetry group As is customary, we label these 
^ G {/, V, T, A, P}, with V and A each grouping together 
four of these one-dimensional irreps, and T six. When 
focusing on a bilinear that transforms under rotations 
as a scalar, vector, tensor, axial vector, or pseudoscalar, 
we shall write for either 1^, 7^, ia^'^ , j^^, or 7|, as 
the case may be. For example, in this compact notation 
the taste-singlet pseudoscalar density is Pj = X7/X- Ap- 
pendix contains explicit formulas for bilinears in the 
taste-singlet representation /, for all F. 

These operators create states such that [45ti47] 



(12) 



where q and q are continuum 4rij-species fermion fields, 
on the right-hand side F is now a (usual) Dirac matrix, 
and ^ is now a four- by- four matrix generator of U(4). 
Together the tensor products ^ T generate U(4n/). 
For nonsinglet ^ €5 T the pseudoscalar meson masses 
depend sensitively on m and a, consistent with chiral 
perturbation theory |Tl]. The fiavor-taste singlet, with 
f (X) T cx l4„j:, should have a mass larger than the rest, 
cf. Fig. JT] but that has not yet been demonstrated nu- 
merically [31] . In this paper, we address this problem by 
studying the eigenvectors of i^stag- 

These lines of theoretical and numerical results lead to 
the picture that 48J 



Dsta,g + m = (^ + 7n) (g) 14 + a^A, 



(13) 



where Ip is the continuum Dirac operator, I4 is the 4x4 
unit matrix, a^A is a taste-symmetry breaking term. 
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and taste-singlet O(a^) corrections are not written out-'^. 
Then 

Det(Atag + m) = [Det(^ + m)]4eT'-'''"'^'^(^+™)"\ 

(14) 

suggesting that 



[Det(£'stag + m)]^/"* ^ Dct(^ + m) x 



i Tr4 lna^A(4»+m)"^ 



(15) 



It is difficult to prove rigorously that the second fac- 
tor becomes benign as a — )■ 0, although a detailed 
renormalization-group argument makes it plausible |5H 
[55]. At nonzero a this factor leads to nonlocality [55] 
(though not the nonlocality discussed in Rcfs. [5U [55]) 
and violations of unitarity. In this paper, we have nothing 
to add to the arguments marshalled elsewhere [TTJ [17^f20] 
that these problems go away as a — 0. 

A separate line of criticism [351 [35] focuses not on 
the ultraviolet taste breaking of gi Tr4ina"A(^-hm)-i 
on the interplay of the rooted determinant with correla- 
tors built from valence propagators. These papers assert, 
without derivation, certain symmetries and properties of 
the 't Hooft vertex that, if true, would imply an unphys- 
ical m dependence of multipoint meson correlators. In 
the next section, we derive, rather than assert, the form 
of the staggered-fermion 't Hooft vertex. Our derivations 
pinpoint where Refs. |351 |36j go astray. Our derivation 
further reveals what is needed for staggered fermions to 
generate the continuum-QCD 't Hooft vertex. Whether 
staggered fermions behave in the needed way depends 
on dynamics, for which a numerical test is needed. The 
(favorable) results of this test are presented in Sec. IV 



III. NEAR-ZERO MODES AND THE 
'T HOOFT VERTEX 

In this section, we discuss the properties of the near- 
zero modes in more detail. We review properties of the 



't Hooft vertex in continuum gauge theory, with one and 
with four species. Then we derive the 't Hooft vertex 
for staggered fermions. We show that the eigenvectors 
must exhibit a certain structure if unrooted staggered 
fermions are to tend to the continuum gauge theory. This 
structure is precisely the criterion presented in Ref. |19j 
for the rooted theory to have a sensible rj' correlator. 



A. Continuum QCD 

In continuum gauge theories, the Dirac operator can 
have genuine zero modes. For a single species, the eigen- 
functions and eigenvalues are denoted ]p(t>a — i^ai'cn 
where A is real, and integer a labels the modes. For 
the modes with nonzero eigenvalue, it is convenient to 
take cr > (cr < 0) for modes with A > (A < 0). 
These modes come in conjugate pairs: = — Ao-, 

= l^<Pa-- III the subspace of zero modes, A = 0, the 
eigenfunctions can be chosen such that 7^(/)t^^ = ±(/)t*\ 
with the integer label t ranging from 1 to k±. For 
n species, the Dirac operator is -^1„, with eigenfunc- 
tions (/)cre^'^\ where the e'^'^^ form an orthonormal basis 
in species space. The number and chirality of zero modes 
is related to the topological charge Q via the index the- 
orem [551 [57] 



(16) 



where n is the number of species, and n± — nk± accounts 
for the species multiplicity. 

The determinant acquires a factor of mass m from each 
zero mode. As m — it would seem that such gauge 
fields would drop out of the ensemble average. But if one 
looks at the eigenvalue-eigenfunction representation of 
the propagator, one finds powers of l/m that cancel the 
powers of m from the determinant. Focusing on |Q| = 1 
and n — 1, so that there is one zero mode, the propagator 
is ((•)|Q|=i denotes average over \Q\ — 1 gauge fields) 



where (f)Q is now used for the zero-mode eigenfunction. One sees that the mode with A = has a canceling factor of 
l/m. The factor (j)Q{x)(j)\{y) is the 't Hooft "vertex" [371[5H]. (If 4>q is localized, as it is around instantons, then the 
"vertex" has support only for x^y near the center of localization.) For the four-point function, there are superficially 



Reference 1481 shows that the Ue symmetry requires the 
dimensions-5 terms of the off-shell Lee-Sharpe [47) effective La- 
grangian to possess coefficients proportional to ma. The appar- 
ent 0(a) effects in Refs. |49l I50| are an artifact of the choice of 
field variables. 
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two powers of 1/m, but two contributions identical apart from their opposite sign cancel each other. This is simply 
the Pauli exclusion principle arising from the Grassmann nature of the fields. 

With n = A fermion species, each mode is replicated four times, so gauge fields with |Q| = 1 yield four zero modes, 
one per species. The determinant yields a factor m*, which is not compensated until the eight-point function: 



with four factors like that in Eq. (17). In higher-point functions, Pauli exclusion again ensures that contributions 
singular in m cancel. Below we are interested in flavor-singlet meson correlators, such as (flavor index contracted; 
Tlf = 4) 

l^Y{^Tfi:{xf)\^2Al^X{{\l + m'f\{4^lTfM^f)\ +0{m), (19) 

\/=l / \f^>0 /=! /|Q|^i 



where the combinatoric factor 24 obtains after cancella- 
tions between many (dis)connected terms. 

Let us now examine a property of the 't Hooft vertex 
that is central to Creutz's arguments [351 [33]. Under the 
anomalous U^(l) transformation 



the n-species 't Hooft vertex transforms as 

n n 

/=1 /=1 



(20) 



(21) 



where the sign is the chirality of the zero mode, j^(f>o = 
±00- If a is a multiple of 27r/n, the prefactor is unity; 
thus, the 't Hooft vertex remains invariant under a Z„ 
subgroup of U^(l) [58] . 

This invariance holds for the full determinant [551 1201 ■ 



Under the transformation (20) with a = 27r/n, one has 



n 

m"("++"-)Det'[^ + me''''^''/"], (22) 

n 

where Det^ denotes the n-species determinant with zero 
modes projected out. The right-hand side follows be- 
cause, by Eq. (16), the phase on the left-hand side is 
trivial. Because e"^/" — e^*27r(n-i)/n^ ^^le twisted mass 

be removed with nonsinglet SU^ (n) trans- 
formations, namely. 



ip 1-^ e' 



V', 



(23) 



where i^ = diag(l, ... 1, — (n — 1)), or any permuta- 
tion thereof. The composition of transformations ( 23 ) 



and (20) with a = 27r/n returns the original determi- 
nant, 77i"("++"-) Det^(.^-|-TO). We have shown here that 
the Z„ in question is not only a subset of the anoma- 
lous U^(l), but also the center of the exact S\J A(n). In 
fact, Z„ is the intersection of the S\J Ain) and Uyi(l)- 



B. Unrooted staggered fermions 

Now we would like to see how staggered fermions repro- 
duce the four-species 't Hooft vertex. Let us now denote 
the eigenvectors and eigenvalues Dstagfs{x) = iXsfsix)- 
We use / for the eigenvectors of £>stag, instead of (j) for the 
eigenfunctions of Ip , because our aim is to study whether 
and how a structure like arises from the fs- As 

before, it is convenient to choose s > (s < 0) for 
As > (As < 0). As mentioned above, the function 
f-s{x) = e{x)fs{x) has eigenvalue A_s = —As, which fol- 
lows from the U^ symmetry. One must bear in mind that 
the relation between eigenvectors f±s originates from a 
different flavor of symmetry than the relation between 
eigen fun ctions (l)±a- In the notation introduced above 



Eq. (12), multiplication by e{x) corresponds to 7p, 



taste nonsinglet that, in a continuum four-species the- 
ory, looks like 7^^^, not 

The first step is to single out the modes analogous to 
the zero modes in the continuum theory. With staggered 
(and most other lattice) fermions, no exact zero modes 
arise, but one expects -Dstag to have some exceptionally 
small eigenvalues [44'. A crisp way to identify them is 
via the spectral flow of the operator [ST] 



stag 



-iD. 



stag +M7/, 



(24) 



with Hsti,gfs{x; fi) = As(/i)/s(x; /i); the eigenvalues of 
-Dstag are iAs(O). From the U^ symmetry, f-s{x; ^i) — 
e{x)fs{x;—^), A_s(/i) = — As(— /x). Near-zero modes are 
those with a nearby zero crossing, A(/io) = for <^ A. 
The (taste-singlet) chirality is then 



Xs = sign A3 (/io), 



(25) 



where the prime denotes differentiation with respect to /i. 
Taking the Ug symmetry into account, we can label the 

positive-chirality modes /|^' with i > ranging from 
1,...,£^ (A,; slightly positive) and i < ranging from 
-~l,...,—£+ (Ai slightly negative). A similar labeling 
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scheme can be adopted for the 
modes /f"''. Note that [61] 



2f_ negative-chiraUty 



(26) 



where Xg is a more common way to identify chirahty ! 62j . 
Modes s and — s have the same value of taste-singlet chi- 
rality (whether defined by Xg or Xs), because 7I implies 
transport over an even number of links and, consequently, 
the e sign factors at the two ends of 7! are the same. 

The spectral flow is elegant but computationally de- 
manding. It is also possible to identify the near-zero 
modes by looking for modes with A sufficiently small and 
X sufficiently close to ±1. Although the spectral flow is 
(presumably) more decisive in borderline cases, in prac- 
tice, especially for the scope of this paper, the compu- 
tational demand seems prohibitive. In Sec. IV we shall 
therefore rely on our experience in Refs. [27, 30j of using 
(A, A") to identify the near-zero modes. 

If staggered fermions generate four species in the 
continuum limit, then the eigenvalues should arrange 
themselves into closely spaced quartets. For nonzero 
modes, four modes should cluster around some distinctly 
nonzero value. For near-zero modes, on the other hand. 



such quartets lie slightly above and below the real axis. 
Ug symmetry dictates that a mode and its e partner have 
the same chirality and, thus, may be assigned to the same 
quartet. If the gauge-field dynamics yield even i± — 2k±, 
then one has quartets. The resulting index theorem is 
then {n± — 4fc±) 



— ri_ = 4(3, 



(27) 



where Q is a pure-gauge definition of topological charge. 
For smooth enough fields and for extensions of Eq. ([s]) 
that smooth out the interaction, both kinds of quartets 
emerge [27l |28] , as does the connection between gauge- 
field topology and the index [27^50] . 

With one fermion field but sets of four near-zero 
modes, the combinatorics underlying the 't Hooft mech- 
anism are less straightforward than in four-species con- 
tinuum theories. Let us focus on \Q\ = 1. Two pairs of 
near-zero modes appear with eigenvalues ±iAi, i = 1,2. 
"Small" means |Ai| ^ (aA)P^A; a power law with p\ = 1 
or 2 suffices, and one expects p\ = 2 |^. Moreover, Xi 
and X2 have the same sign (with several actions pT l lSO ] ). 
and we shall see in Sec. IIVI that these features also 
hold for the highly-improved staggered-quark (HISQ) ac- 
tion IMl. 



To derive the 't Hooft vertex explicitly, let us examine the (fermion) eight-point function, which for staggered 
fermions is 



\/=i 



\i=l 



s>0 



(28) 



'|Q|=1 



where the propagator 

Gix, y) = {x{x)x{y))x,x - E ^-"^T^ (29) 



all 



with the sum running over near-zero and nonzero modes. 
Neglecting in Eq. (28) the near-zero Ai relative to m, the 
near-zero-mode terms contribute to Eq. ([28|) as 



m'^ det G{xf,yg) 



det Mxf)det f}{yg) + 0{m), (30) 



where i,j G {—2, —1, 1, 2}. In higher-point functions, the 
Pauli exclusion again ensures that contributions singular 
in m cancel. 

The product of determinants on the right-hand side 
of Eq. (30) is the 't Hooft vertex for (unrooted) stag- 
gered fermions. To repr oduce the product of four fac- 



tors of (po(l>Q in Eq. (|18|), the four staggered eigenfunc- 



tions fi, i e {—2,-1,1,2}, must have structure sim- 
ilar to 4>o{x)e''^\ One could seek such structure in a 
basis where a taste index looks obvious, but because 
taste is, fundamentally, a quantum number of the shifts. 



i.e., single-link translations, gauge-dependent roughness 
of the gauge field would obscure it. 

The way forward is to contract the x X fields 
into color singlets. The contractions must also be taste 

singlets, because a nonsinglet corresponds to e'^'^^^e'--'-'. 



^ ^ I4, which need not vanish when j ^ i. In Eq. ( |28[ ) we 
thus replace x{^f)x{yf) with a taste singlet xXjx{xf)- 
Contracting Eq. ( 30 ) in this way, one is led to consider 



flTif,{x) 



(31) 



with, recall, some parallel transport implied by Tj. The 
't Hooft vertex simplifies in the desired way if 



C^.(x)cx<5,,[l + 0(aPcr)], 



(32) 



for Tj = Ij, '-fj, iaj . If, further, the proportionality 
fulfilled by an (approximately) i-independent diagonal 

Cii{x), Cfj would then mimic ^jFc/ioe^'^ e^^^ oc Sij. Ap- 
proaching this limit as a power law with p ^r = 1 or 2 

presents 



IV 



suflices, and one expects p^r = 2 [48]. Section J 
numerical results for these local overlaps, including a de- 
pendence. 
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For Tj = 7/ , 7j^, the local overlaps (fj behave some- 
what differently. In continuum gauge theory, the zero 



modes satisfy (pi"^^ ' 7^(/)i 
7^ anticommutcs with 7^ and 7 
spin and taste degrees of freedom emerge from staggered 
fermions via the same dynamical mechanism, so the di- 
agonal and Q: should vanish commensurately with 



T 



= 0, because 



The 



the off-diagonal Cj,, F/ 



1/, ll 



The local overlaps of continuum nonzero modes of dif- 
ferent species also vanish (trivially, because e'^'^^-'^e^'^^^ = 
5riT2)- Therefore, within a quartet of staggered- fermion 
nonzero modes, continuum QCD is reproduced if Cr^, 
r ^ also vanish as a 0. 



Assuming Eq. (32 1 holds, it is easy to see that 



\{x^ifX{xf)) = ( XliXl + m-'f ^ui^^)(:]!i^2)Cll{x^)<:l^{x,)\ +0(m), (33) 

\f = l I \<7>0 (tjkl) 



where the sum runs over the 4! — 24 ways of choosing 
distinct (ijkl) from {—2, —1, 1, 2}. 

Let us now discuss the Z„ (now Z4„^ ) symmetry men- 



tioned at the end of Sec. Ill A The anomalous 11^(1 



and most of the softly broken, nonanomalous 811^(471/ ) 
emerge only in the continuum limit. Some passages in 
Refs. [35l[36] seem to assign a pertinent role to the Ue(n/) 
symmetries, which are exact even at nonzero a. These 
symmetries are a distraction at best: the group Ue(ny) 
intersects with the relevant Z4„^ , which is the center of 
SU(4n/), only at —l^nf 



C. Rooted staggered sea 

With rooted staggered fermions, two changes are car- 
ried out. In addition to using the rooted determinant ([2]), 
the simple combinatorics of det^j^g) G{xf,yg) must also 
change |39]. For example, the taste-singlet pseudoscalar 
meson propagator is replaced with 

{X7!x{x)xih{y))u ~lC{x,y) + ^,D{x,y), (34) 

where the connected and disconnected contributions are 

C(x, y) = (D tr [7! G(x, yh'jGiy, x)])^, (35) 
D{x,y) = (Dtr[7fG(^,a:)]tr[7fG(y,2/)])^, (36) 

where D is the rooted determinant ([2|, the trace is over 
color, and the translations implied by 7I act to the right 
(left) on the first (second) argument of G. The correlator 
in Eq. (34) couples to the analog of the flavor-singlet 77' 



meson in QCD, and similar constructions hold for other 
taste-singlet bilinears. 

The combinatoric factors in E q. ([M] ) follow immedi- 
ately from considering [HI (THl 1301 Si] 



Det [{p + to) 



n/4 



(37) 



where — inside the braces — one has four copies of non- 
controversial fermions. Equation (371 together with a 



source for a single species provide an engine to generate 
the combinatorics of rooting (in general): to obtain n 
species from 4, a term with t traces over color receives a 
factor 53 



(-i)' 



(38) 



For Eq. (37 1 to be relevant to staggered fermions. 



the dynamics must ensure Eq. ( 13 ) and, in particular, 
Eq. ( 32 ) , as we now show. The single-flavor determinant 
becomes (for \Q\ = 1) 



(39) 



s>0 



Neglecting compared to to again, the first product 
collapses to \m\. The near-zero- mode contributions are 
then 



C{x,y) 



D{x,y) - 



2J 



(40) 



|Q|=1 



<l{^)q,{y)) , (41) 

/|Q|=1 



where D' is the s > product in Eq. (39), and i,i £ 



{-2, -1, 1, 2}. If Eq. ([321 holds, then the sum in Eq. ^ 
collapses to terms with i = j, apart from lattice artifacts. 
Thus, C has 4 contributions singular in 1/|to|, whereas 
D has 16. With the correct combinatoric factors, they 
cancel. 

It is, perhaps, instructive to exhibit the three-point 
correlator. Assuming Eq. ( 32 1 and homing in on the zero- 
mode contributions, 
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X^nxixi)xri2X{x2)xri3Xix3) - 3{tr[r/iG(xi,a;2)r72G(a;2,X3)r/3G'(x3,xi)] + 1 perm} 

+ ^{tr[r/iG(xi,xi)]tr[r/2G(x2,a;3)r/3G(x3,a;2)] +2 perms} 

- ^ tr[rnG(a;i, xi)] tr[r/2G(a;2, X2) tr[r,3G(a;3, 2^3)] (42) 

^ M (_2 + 3 - 1) {xi)e^ (X2)C^^ {X3)C^^ {x,), (43) 



where |m| comes from the rooted determinant. Here sums 
over the four staggered-fermion near-zero modes cancel 
the exphcit factors of \. The |TO|/m^ contributions can- 
cel in a similar way. Earlier work [18l |40l [41], tacitly 
assumed Eq. (32 1; in particular, Ref. [H] shows how the 
combinatorics work for higher-point 't Hooft-vertex ef- 
fects. 

In Refs. |35J UHl El], Creutz disregards the cancella- 
tions stemming from the correct weighting of different 
contributions to flavor-taste-singlet correlators. He con- 
siders more primitive combinations, like any individual 
line in Eq. (42 ), which clearly are singular as to — 0. He 
then draws two incorrect inferences. First, he claims that 
the normal cancellations connected with Pauli statistics 
cannot arise. Combining the correct weights with the 
assumption (tested below) Eq. ( 32 1 , one sees that this is 
not the case. The outcome is not too mysterious: as taste 
emerges into a species-like quantum number, the correct 
set of correlators averages over them. 

The other misstep is to assert that the Z4„j: symme- 
try of the unrooted 't Hooft vertex cannot be reduced 
to Z„j,. This is incorrect, because, while the rooted de- 
terminant clearly retains the symmetries of the unrooted 
determinant, the 't Hooft vertex stems from the com- 
bined behavior of determinant and valence propagators. 
The replacement of the combinatoric factors of traces 
with (38) effectively projects the symmetry emerging in 
the chiral limit from SU(4n/) [taking n 



(38)] 



to 
the 



= 1 in 

SU(n^). Since the relevant symmetry is the center o 
emergent flavor symmetry, one has . 

Many of these points have been made before [TSl 1301 
I4TI [65] , but until now it has always been assumed that 

(Refer- 



(32) 



the tastes decouple as posited in Eq. 
ence (TH] noted the necessity of this assumption.) Our ap- 
proach can easily be extended to taste-nonsinglet flavor- 
singlet correlators, and the properties of the local over- 
laps with nonsinglet Fs will not enjoy the cancellation. 
We shall now compute the Cfj nonperturbatively, to flnd 
out whether the tastes couple to each other at the strong 
scale Aqcd or at the cutoff scale a~^. 



IV. NUMERICAL RESULTS 

In this section we present our numerical methods and 
results. First we explain the motivation for studying 
improved discretizations and why it suffices to compute 



TABLE I; Details of the gauge configurations used: /3 is the 
bare gauge coupling, a the lattice spacing [68], V the space- 
time volume in lattice units, and L the linear size in physical 
units. The final column gives the number of configurations 
in each ensemble with \Q\ = 1. We refer to set 1 as having 
a "coarse" lattice spacing, sets 2, 3, and 4 as "intermediate", 
and set 5 as "fine". 



Ensemble 




a (fm) 


V 


L (fm) 


#{|QI = 1} 


1 


4.6 


0.125 


12* 


1.50 


294 


2 


4.8 


0.093 


12* 


1.12 


806 


3 


4.8 


0.093 


16* 


1.49 


424 


4 


4.8 


0.093 


20* 


1.86 


288 


5 


5.0 


0.077 


20* 


1.54 


430 



their eigenvalue spectrum on quenched gauge fields. We 
present results for eigenvalues and chirality with the 
HISQ action. These results are qualitatively similar to 
those obtained with the Asqtad and Fat7x Asqtad actions 
in Refs. [271 [3D]) so we focus here on |Q| = 1. Then we 



show results for the overlaps, Q, defined in Eq. (31), and 
test their behavior as a function of lattice spacing against 
Eq. ( 32 ) . Finally we discuss correlators for mesons of 
different in turn, starting with pseudoscalars where 
the issues are particularly important. Taken together, 
our results demonstrate how the behavior of the different 
contributions from near-zero and nonzero modes matches 
that expected in the continuum. 



A. Methods 

In this paper, we use the same ensembles of SU(3) 
gauge fields as in earlier studies of eigenvalues and chi- 
rality [271 [SO]- They are quenched configurations, omit- 
ting the effects of sea quarks. They are generated with 
a Symanzik-improved gauge action, so that the tree-level 
errors are removed [66] . and tadpole- improved cou- 
plings in this action, so that loop corrections are re- 
duced 67^ . Three different values of the gauge coupling 
are used, giving three widely separated values of the lat- 
tice spacing, covering the range of typical unquenched 
lattice-QCD calculations [TT], so our results should per- 
tain directly to them. At the middle value of the three 
lattice spacings, we have three different-sized lattices in 
order to check the volume dependence. The parameters 
for the configurations are given in Table [l] 
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It is sufficient to study these issues in the quenched 
approximation, because we aim to test a structural prop- 
erty of staggered fermions in fixed-Q sectors. In par- 
ticular, omitting the determinant decouples Creutz's in- 
frared concerns from others' ultraviolet concern that 
taste breaking remains in the continuum limit. If the 
eigenvectors satisfy Eq. (32) strongly enough, then the 



't Hooft vertex and the consequent cancellation of mass- 
singular contributions to the connected and disconnected 
flavor-singlet meson correlators should work out in gen- 
eral. We shall see that this is the case. 

With the original staggered- fermion action, Eq. (|3]), 
the interaction connects adjacent sites. Very large dis- 
cretization errors arise in a wide range of observables, 
washing out the expected quartet structure in the eigen- 
value spectrum. These discretization errors have been 
traced to taste-changing interactions from gluons with 
one or more components of momentum 7r/a 
Because of the gluon exchange, these effects are formally 
of order a^a^, i.e., as times smaller than normal dis- 
cretization effects [7D]. In order to reduce these taste- 
changing effects, it is necessaryto smear the gauge field, 
replacing J/^ and Uj^ in Eq. i3« with sums of products 
of link matrices tracing out more complicated paths be- 
tween X and X ± fia [521 US IZD ■ 

Several staggered-fermion actions have been developed 
along these hues. The Asqtad [73] and Fat7x Asqtad [30] 
actions exhibit a reduction, relative to the nearest- 
neighbor action in Eq. (|3|, in splittings between pseu- 
doscalar mesons of different taste [26,i30j. Similarly, with 
these actions the quartet structure of the eigenvalue spec- 
trum more clearly emerges [27l [30] . 

Here we have calculated low-lying eigenvalues and 
eigenvectors for the highly improved staggered-quark 
(HISQ) action [63 , reusing the same gauge-field config- 
urations. The HISQ action supersedes the FatTx Asqtad 
action; it is essentially the same but corrects the smear- 
ing at the second stage to remove fully the discretiza- 
tion errors that the smearing introduces. As we shall 
see in Sec. |IVB| this change makes only a small effect. 
The eigenvalue quartet structure is very clear with the 
HISQ action, which is reflected in other properties that, 
by now, have been thoroughly tested: small pseudoscalar 
mass splittings and small discretization errors, even for 
heavy quarks [Ml [UlIS] ■ 

Appendix [B] provides explicit equations for the 
smeared actions. 

We use the Lanczos algorithm to calculate the low- 
lying eigenvalues, iX, of the anti-Hermitian massless 
HISQ Dirac operator, I^hisQi defined implicitly in 



Eq. (B6). Owing to its red-black checkerboard structure. 



the calculations can be simplified by using the Hermi- 
tian positive semi-definite operator — ^hisqi projected 
onto either the red (even) or black (odd) sites of the 
lattice. This yields A^, from the smallest values up- 
wards, and eigenvector /, on the chosen half of the lat- 
tice. The eigenvalues of Dhisq are then zLiX, and the 
corresponding eigenvector on the other half of the lattice 



is ±_Dhisq//*A. This construction automatically imple- 
ments the requirement that the eigenvectors correspond- 
ing to eigenvalues iX and —iX are simply related by mul- 
tiplication with e{x). Thus, on the odd (even) sites, the 
— sth eigenvector is opposite (same) in sign as the -|-sth 
eigenvector. 



B. Eigenvalues and chirality with HISQ 

Figure [2] shows the four near-zero eigenvalues as well as 
the 16 pairs of nonzero eigenvalues of -Dhisq with small- 
est |A|, obtained on typical \Q\ = 1 configurations from 
ensembles labeled 1 (coarse), 3 (intermediate), and 5 
(fine) in Table |l] These lattices have similar physical vol- 
ume but lattice spacing varying from 0.125 to 0.077 fm. 
The anticipated picture is unmistakable: four (and only 
four) very small eigenvalues appear, followed by distinct 
quartets. As the lattice spacing decreases, eigenvalues 
within a quartet come closer and closer to being degen- 
erate, typically by forming two close-by almost degener- 
ate pairs. The near-zero modes are typically, on these 
lattices, at least an order of magnitude smaller than the 
low-lying nonzero modes. 

The Lanczos algorithm also gives the eigenvectors cor- 
responding to these eigenvalues. Normalizing them to 
have modulus 1, we compute the chirality X in Eq. (26 1, 



using the smeared Wfj, matrices [Eq. (B5)] instead of f/^. 
Reference [3D] showed that it makes little qualitative dif- 
ference to the results whether the original C/^;, Asqtad 
Vf, [Eq. (|B2])], or Fat7x Asqtad [Eq. ([bs])] are used. 
The numerical values of the chirality may change, but 
the picture remains qualitatively the same. 

Because lattice artifacts break the taste-singlet sym- 



0.25 
0.20 
0.15 
0.10 
0.05 



Xa 0.00 

-0.05 
-0.10 
-0.15 
-0.20 
-0.25 



1 1 1 1 1 T 

• a = 0.125fm, y= 12 

■ a = 0.093fm, ^=16* 

* a = 0.077 fm, ^=20* 



"1 — r 



"1 — n 



J L 



J L 



J L 



1 2 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
lil \s\ 

FIG. 2: The four near-zero eigenvalues (left panel) and the 
16 lowest-lying nonzero pairs of Dhisq eigenvalues on a typ- 
ical |Q|=1 configuration from sets 1 (red circles), 3 (green 
squares), and 5 (blue triangles). For clarity, some modes are 
offset horizontally. 
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metry, the chirality defined in Eq. ( 26 ) takes values that 
are not simply 1 



and [62]. References [27l [30] found, 
however, that it is easy, especially with improved gauge 
and staggered-fermion actions, to separate the near-zero 
modes with relatively large chirality, close to 1, from the 
other modes with chirality close to 0. The number of 
near-zero modes defined this way agrees with the index 
theorem, Eq. ( 27 ) , and pure-gauge definitions of the topo- 
logical charge. 



The agreement between the index and 
the gauge-field topological charge improves as the lattice 
spacing gets smaller. On the a = 0.077 fm ensemble, 
the disagreement for Asqtad and Fat7xAsqtad is just 
2% [3D], which is no worse than the ambiguity between 
different gluonic definitions. For this paper, we therefore 
simply take the index to classify the topology. 

Figure [3| shows the chirality values for the HISQ ac- 
tion versus eigenvalues on all configurations defined to be 
of topological charge ±1 via the index. To reduce clut- 
ter. Fig. [3] shows only the two near-zero modes and the 
two lowest-lying nonzero modes. (Because X^g = Xg, 
we count only the positive-A modes here.) One sees a 
clear separation of large and small chirality values, es- 
pecially so on the finer configurations. Although the 
values corresponding to the maximum chirality do not 
change very markedly from coarse to finer lattices, the 
spread of results becomes much narrower. The small chi- 
rality values, corresponding to nonzero eigenmodes, fall 
rapidly to zero with lattice spacing. We take X > 0.4 
(drawn on the graphs) to indicate large chirality and then 
count the number of eigenvalues (with positive A) that 
have large chirality. Configurations with two (positive-A) 
large-chirality modes are taken to be \Q\ = 1 configura- 
tions. Table |T| lists the number of such configurations for 
each ensemble. More general scatter plots with results at 
IQI > 1 and the Asqtad and Fat7x Asqtad actions have 
been given in Ref. [30j . and with HISQ look very similar. 



\X\ 



a = 0.093 fm 



a = 0.125 fm 



Xa 



FIG. 3: The absolute value of the chirality \X\ plotted against 
eigenvalue, Xa, in lattice units for the four lowest (positive) 
eigenvalues for the |(5| = 1 configurations in ensembles 1, 
3, and 5. The dotted line on each graph indicates \X\ = 
0.4, which is used to separate large and small chirality in 
determining the value of Q (see text). 



Results for 



Using the eigenvectors determined in the previous sec- 
tion we now go on to look in more detail at the overlaps 
of the near-zero-mode eigenvectors that are relevant to 
the 't Hooft vertex. Figures |4|-[6| show scatter plots and 
histograms of the (^fj distributions for F = 1, 7^, and 7'^, 
and i,j = ±1,±2 on one or two configurations, ranging 
over all x. Each figure displays this information, from 
top to bottom, for the coarse (a — 0.125 fm), intermedi- 
ate (a = 0.093 fm), and fine (a = 0.077 fm) lattices, at 
(nearly) fixed physical volume (sets 1, 3, and 5). The four 
panels in each case show the scatter of (fj in the com- 
plex plane (upper left), the histogram for ReCfj (lower 
left), the histogram for Im^'^- (upper right), and the his- 
togram for (lower right). The number of points in 
the histograms for each set is the lattice volume, V, of Ta- 
ble [l| Note the logarithmic scale on the histogram plots. 
Red points and lines denote diagonal C,fi, and black off- 
diagonal (fj (j ^ i). In the case of the vector overlap, 
we separate the off-diagonal (fj into two. Black 



=/= \i\ and the case of j — —i is shown 



F = 7' 

is reserved for 
in blue. 

The most striking feature for the scalar (Fig. |4]) and 
pseudoscalar (Fig. [s]) is how different the diagonal and 
off-diagonal distributions are. The diagonal scalar over- 
lap is a sum of absolute squares, so it is real and pos- 
itive. Because f^i{x) — e{x)fi{x) and the taste-singlet 
scalar operator is local, is equal to Qj^ on even sites 
but real and negativ e on odd sites. Upon averaging over 
a hypercube in Eq. (Al), cancellations render Q _^ rela- 
tively small. It is visible on Fig. |4|as a black line stretch- 
ing along the negative real axis; the positive part being 
invisible underneath the red line for Q.^. The off-diagonal 
(bl 7^ Nl) scalar overlap is a complex number of ran- 
dom phase. The width of all the histograms falls going 
down the column of plots as the lattices become finer. 
What is crucial for the taste structure of the eigenvec- 
tors, however, is the relative width of the histograms for 
ICj'jl for i ^ j compared to that for From the plots 
it can be seen that the width of the off-diagonal distri- 
bution is falling faster with lattice spacing than that of 
the diagonal. Figure |4| shows a single configuration with 
IQI = 1, but we have examined others, and they look the 
same. 

Figure [sj for the pseudoscalar case shows two config- 
urations, one each with Q = -1-1 and —1. The plots 
behave in the same way as the scalar overlaps, except 

that is real and negative for Q = —1, as a conse- 
quence of parity. From the same arguments as above, 
since the taste-singlet pseudoscalar operators links odd 

sites to odd sites and even sites to even sites, C7-i is also 

'5 

real and takes the same or opposite sign to on odd 
sites or even sites. Q _i is therefore not visible beneath 



on Fig. [sj Once again, looking down the plots, we 
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FIG. 4: d^lj on coarse (top), intermediate (middle), and fine 
(bottom) IQI = 1 gluon field configurations, with j — i (red) 
and j i (black), i,j — ±1,±2. Note the logarithmic j/-axis 
scale for the histograms. 



FIG. 5: (^^ on coarse (top), intermediate (middle), and fine 
(bottom) IQI = 1 gluon field configurations, with j — i (red) 
and j ^ i (black), i,j = ±1,±2. Note the logarithmic y-axis 
scale for the histograms. 



12 



Re(g 

-0.0002 0.0002 



100 10000 







-0.0002 0.0002 

Re(Q 



Re(i;) 

-0.0002 0.0002 



-0.0002 0.0002 






-0.0002 0.0002 

Ee(i;) 



ReCC) 

-0.0002 0.0002 



-0.0002 0.0002 



100 10000 





-0.0002 0.0002 

Re(i;) 




FIG. 6: (^7j on coarse (top), intermediate (middle), and fine 
(bottom) [Qj = 1 gluon field configurations, with j — i (red), 
j = -i (blue), and \ j\ ^ \i\ (black), i,j = ±1,±2. Note the 
logarithmic y-axis scale for the histograms. 



see clearly that the width of the off-diagonal distribution 
(combining j = —i and |j| 7^ \i\) decreases with lattice 
spacing, relative to the diagonal distribution. 

With (^j the behavior differs. Recall that should 
vanish for all even j = i. As seen in Fig. [m we find 
Q^ to be pure imaginary, which follows from the defini- 
tion of the operator 7}', Eq. ( A2 ); we find (^^__^ to be pure 



real, which follows by changingthe sign of the odd pieces 
of f-i{x) relative to fi{x), because the vector operator 
couples even to odd sites and vice versa; and we find (^j , 
\j\ 7^ \i\, to be complex and of random phase. In this 
case, however, the widths of all three distributions not 
only are the same (when nonzero) but also decrease with 
decreasing lattice spacing together. Indeed the widths of 

all I distributions are similar to the widths of the 

5 

and iCi^- I distributions, j 7^ i. 

To visualize the lattice-spacing dependence more di- 
rectly, we plot in Fig. [7] the width of the (fj distribu- 
tions, appropriately normalized, vs a^. The widths are 
defined by the central 66% of the data in the lower right 
histogram for \Cfj\, but calculating this histogram for ten 
configurations instead of just one or two. The errors are 
estimated by comparing the widths for two subsets of 
five configurations. The values we obtain for the widths, 
and their errors, are given in Table |Tl] Since the eigen- 
vectors are normalized to have modulus 1 at each lattice 
spacing, the widths do not have a physical interpretation. 
The best that one can do is to normalize the off-diagonal 
widths against diagonal widths, as is appropriate for the 
test of Eq. (32). This ratio of widths is plotted for the 
scalar and pseudoscalar in Fig. [7| For the vector, we 
have no diagonal quantity that survives in the continuum 
limit, so we normalize instead against the diagonal pseu- 
doscalar width. Although it is difficult to be quantitative 
(full ensemble averages of the widths are too costly, and 
the determination of the lattice spacing in the quenched 
approximation is ambiguous) , the trend in Fig. [t] is clear 
and consistent with what is needed according to Eq. ( 32 1 



Figure [7] shows, with dashed lines, representative fits 
as a polynomial in to our results. The fits include a 
constant plus quadratic, quartic, and sixth powers of a. 
The slope of the nth polynomial term is constrained by 
a Bayesian prior to a size of (1.0 GeV)" suggested by the 
slope of pion taste splittings [7S] . It is very easy to obtain 
good fits with any combination of different polynomials, 
for example, including or not including a linear term, so 
it is not possible to say definitively what the lowest power 
of a is that appears in the a dependence of the Qj . The 
solid points on the plots in Fig. [7] give the a = value 
of the width ratios, compatible with zero in all cases. 
Thus, our results are consistent with the expectation in 
Eq. (|32| , although the data are not able to determine p^r 
in a definitive way. 

We have also investigated the volume dependence of 
the (^fj for i,j ~ ±1,±2 at the intermediate lattice spac- 
ing (i.e., on sets 2, 3, and 4), and these results are also 
included in Table In) We see that the widths again fall 
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TABLE II: Widths, in lattice units and multiplied by 10*, obtained for different \(^fj\ histograms on each set of \Q\ = 1 gauge 
configurations. The first column gives the set and then subsequent columns list the width, with an estimate of the error, for 
different F and i,j combinations, i,j = ±1, ±2. 



Set 



7 



T 



-J r 



^ lil 1, 



1, i/i 



1 3.93(31) 0.296(26) 0.242(11) 0.202(21) 0.224(1) 

2 4.08(56) 0.180(11) 0.188(23) 0.146(12) 0.155(16) 

3 1.49(3) 0.0635(38) 0.0567(36) 0.0541(42) 0.0561(22) 

4 0.548(26) 0.0390(16) 0.0242(10) 0.0246(17) 0.0258(4) 

5 0.716(14) 0.0230(7) 0.0211(2) 0.0206(2) 0.0219(5) 



6.84(17) 0.285(23) 
6.07(70) 0.158(9) 
2.15(6) 0.0566(27) 
0.846(29) 0.0406(22) 
0.954(16) 0.0199(5) 
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7: (color online) The top plot shows the width of \C^. 
i, divided by the width of \Cfi\, for F = 1 



±1, ±2. The F = 1 case is given by open circles (red online), 
along with a representative polynomial fit in the square of the 
lattice spacing and the corresponding value in the continuum 
limit (filled circles). The equivalent results for F = 7'' are 
given by open and closed squares (blue online). The lower 

plot shows the width of | divided by the width of |("^ | 
plotted against the lattice spacing. The case i — j is given 
by (black) open and closed squares, the case i = —j by (red) 
open and closed circles, and the case |i| 7^ |j| by (blue) open 
and closed triangles. 



as the volume of the lattice increases. Naively this is 
simply a result of the normalization of the eigenvectors 
to 1 over an increasing number of lattice sites. Indeed 
the widths do seem to have simple behavior, inversely 
proportional to 1/V, at least for the diagonal scalar and 
pseudoscalar widths and the vector widths. Note that 
this is not inconsistent with the fact that, for example, 
the very high values of the pseudoscalar diagonal over- 
laps are localized around the instantons that give rise to 
the near-zero modes. 

The pseudoscalar and scalar widths behave quite dif- 
ferently as a function of lattice spacing than they do as 
a function of volume. We can see this by comparing the 
histograms in Fig. [s] for (f^ on the fine lattices, set 5, 
and the large intermediate volume lattices, set 4. Both 
of these have 20'* lattice points. We see that the diagonal 
distribution is broader on the finer lattices and the off- 
diagonal distribution markedly narrower, consistent with 
the fairly rapid fall with lattice spacing of the ratio of 
the widths seen in Fig. [7] For the vector case, as is clear 
from Table [nj the behavior of the widths with the lattice 
spacing is only slightly steeper than that with volume. 
However, this still represents a fall to zero with lattice 
spacing when compared to the diagonal scalar and pseu- 
doscalar overlaps which survive the continuum limit, as 
we see in Fig. [7) 

We have not shown histograms for the axial vector or 
tensor operators. We have looked at these operators in 
terms of the relevant meson correlators (see the following 
subsections) and they give qualitatively identical results 
to the vector and scalar/pseudoscalar cases, respectively. 
It therefore seems unlikely that they would upset the 
picture gleaned here. 



D. Flavor-singlet meson correlators 



Our results in Sec. |IV C| show how the taste-singlet 
overlaps, C,fj{x), of different near-zero-niode eigenvectors 
behave as expected to give the correct continuum be- 
havior for the 't Hooft vertex. Here we show explicitly 
how this translates into the correct continuum behavior 
for the near-zero-mode contribution to the fiavor-singlet 
meson correlator. We also look at nonzero-mode contri- 
butions, as well as flavor-nonsinglet correlators, wherever 
they are useful to fill out the picture obtained. As dis- 
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cussed in Sec. Ill it is sufficient to work in the quenched 



approximation since the structural issue of the behavior 
of the eigenvector overlaps — in a fixed- sector — is the 
same whether sea quarks are included or not. 

To relate results as closely as possible to those of a 
complete meson correlator calculation in lattice QCD, 
we consider meson correlators projected onto zero spatial 
momentum by summing over spatial sites. This leads us 
to consider a modification of the eigenvector overlaps 



(44) 



summing over a time slice instead of a 2** hypercube. 
Then the zero-momentum connected and disconnected 



contributions can be constructed, as in Eqs. (40) and (41 ) 
from correlations of time-slice overlaps 



{t)clit + T), 



-T). 



(45) 
(46) 



well as near-zero modes. 

The full connected correlator is 
j:^^yC{x,t;y,t + T), with C{x,y) 
Eq. (35). Similarly, the disconnected 
D{T) = D{x, t;y,t + T), following Eq. p]). C(r) 

is then 



C{T) ^ 
defined in 
correlator is 



made 



up of Xrs correlated overlaps (on the 
quenched configurations that we are studying) as 



C{T) 



Xrs{T) 



{iXr + m){iXs + to) 



(47) 



where we have made explicit the dependence on the 
eigenvalues in the denominator. For D{T) we have 



DiT) 



YrsiT) 



{iXr + m){iXs + to) 



(48) 



Note that X^j.{T) = Y^^{T) by construction, and we con- 
sider values for r, s that correspond to nonzero modes as 
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FIG. 8: Histogram of 1^^^ j, for i — j (top) and i j (bottom), 
comparing results on the fine lattice (set 5 — red/gray) and the 
large intermediate lattice (set 4 — black). 



The disconnected correlator factorizes into the product of 

— r 

sums over diagonal overlaps C^^, but the connected cor- 
relator contains overlaps between different eigenvectors. 

Note that the factor \m\. 
the Hf 



of Eqs. (|35| and (|36|) from 
af- 



1 sea quark determinant is missing. T 
fects the weighting of the particular configurations in the 
ensemble and therefore the quantitative results obtained 
for C{T) and D{T). However, it does not affect qualita- 
tively the properties of the Xrs factors that we demon- 
strate here, which are evident in a fixed-jQI sector and 
even, in some cases, on a configuration-by-configuration 
basis in their contribution to C(T) and I?(T). 
As discussed in Sec 



III 



we then have to test whether 
the near-zero modes give rise to a divergence in the cor- 
relator for flavor-singlet meson 7? as to — >■ 0, when the 
connected and disconnected contributions are combined 



with their appropriate taste factors of 4 and 16 [Eq. (34) 



M"{T) = (X^(r))^ = \C"{T) Y^D"{T). (49) 

To obtain a finite result as to — for M{T) we need 
the near-zero-mode contributions to cancel between C (T) 
and D(T). This in turn requires the off-diagonal corre- 
lated overlaps, Xij, i ^ j, between different near-zero 
eigenvectors in the same staggered eigenvalue quartet to 
vanish in the continuum limit. Then each quartet be- 
haves as four copies of a single mode and, including the 
factors of 1/4 and 1/16, reproduces within M{T) the be- 
havior expected of eigenmodes of the Dirac operator in 
the continuum. We show how this works explicitly for 
the examples of scalar, pseudoscalar, (axial) vector, and 
tensor mesons in the following subsections. We do this 
with the same |(5| = 1 quenched configurations used in 
the previous subsection. The correlator results are, how- 
ever, averaged over all \Q\ — 1 configurations for each 
ensemble, rather than just 10. 
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FIG. 9: Diagonal and off-diagonal correlated overlaps (T) 
between eigenvectors 1 and 3 that contribute to the discon- 
nected piece of the 77' correlator at zero spatial momentum. 
Results are given for the average over |Q| = 1 configurations 
in the fine ensemble, set 5. The lower plot shows the off- 
diagonal correlated overlap summed over time separation as 
a function of the square of the lattice spacing. 



FIG. 10; The correlated overlaps X^- (T) between near-zero 
modes 1, 2, and —1 that contribute to the connected piece 
of the Tj' correlator at zero spatial momentum. Results are 
given for the average over |Q| = 1 configurations in the fine 
ensemble, set 5 (top) and the coarse ensemble, set 1 (bottom). 



1. Flavor-singlet pseudoscalar mesons 

We first discuss the important case of the calculation 
of the flavor-singlet pseudoscalar meson (77') correlator, 
and the associated case of the flavor-nonsinglet meson 
(tt) correlator. In continuum QCD with, say, two equal 
mass light quarks, this is readily analyzed in terms of 
the eigenvectors and eigenvalues of the massless Dirac 
matrix. For |(5| = 1 gluon field configurations there is 
one zero mode with a chirality (j)Qj5(f>o ~ ±1- The tt 
meson correlator has no disconnected contribution and 
the connected contribution is readily seen to obey, on a 
given gluon field configuration. 



T 



1 



(50) 



where m is the quark mass and the sum is over all eigen- 
modes of the massless Dirac matrix, including the zero 
mode. Each eigenmode contributes a correlated overlap 



of 1 when summed over T. For the zero mode, this comes 
from the square of the chirality. The nonzero modes have 
chirality zero but still contribute a correlated overlap fac- 
tor of 1 because the 75 matrix connects the modes a and 

— cr with A_CT — —Act. Then (continuum) C^. _g.{t) = 1 
from eigenfunction normalization. 

The 77' meson correlator is made from the same eigen- 
vectors but now has a disconnected contribution com- 
ing from the zero mode that exactly cancels the zero- 
mode contribution to the connected correlator. Thus, on 
a given configuration. 



1 



(51) 



There is now no contribution from the zero mode and the 
correlator is finite as m — ?> 0. This continues to be true 
on averaging over gauge fields and including determinant 
factors. 

Now let us show how staggered fermions reproduce 
Eqs. (501 and (51 1. For the Goldstone vr meson cor- 
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FIG. 11: The correlated overlaps X^^ (T) between nonzero 
modes 3 and 4 that contribute to the connected piece of the 
t;' correlator at zero spatial momentum. Results are given 
for the average over |Q| = 1 configurations in the fine ensem- 
ble, set 5. Results are also shown for the correlated overlap 
between near-zero mode 1 and nonzero mode 3. 



3.- 



X 3.-6 




FIG. 12: The correlated overlaps X^a {T) between nonzero 
mode 3 and modes —3, —4, —5, and —6 from its Ue mirror 
quartet. Results are averaged over \Q\ — 1 configurations for 
set 5. 



relator, it is straightforward and mechanical 
r = 7p = e(a;) connects eigenvectors fs and /_s 
correlated overlap contribution is again 1, when 
over T, simply from eigenvector normaUzation. 
ference with the continuum case is that this is 
for the near-zero modes. Thus, we obtain an 
very similar to that in the continuum on a sin; 
field configuration: 



Then 
and the 
summed 
The dif- 
also true 
equation 
gle gluon 



E 



+ m 

1 



q 1 



+ 0(a2). 



(52) 



where s is a sum over all modes including the near-zero 
modes and the factor of 1/4 is the same as in Eqs. (34) 
and (|49|). Since the staggered eigenvalues come in quar- 



tets that become degenerate in the a — ;> limit, the lower 
equation replaces the 4 eigenvalues in a quartet by their 
mean square and sums all quartets, q, including the near- 
zero- mode quartet. This then clearly reproduces the con- 
tinuum Eq. ( 50 ) as a — )■ 0. 

The flavor-singlet correlator is constructed differently 
and includes both connected and disconnected contribu- 
tions. With staggered fermions, we must use the flavor- 
taste-singlet pseudoscalar, T = jj. In demonstrating 
that Eq. (51) is reproduced, we also show that the corre- 



lated overlaps behave so as to give a finite result for the 
77' correlator. 

From our earlier results on chirality, we can anticipate 
what the disconnected correlated overlaps look like. 



Because J2t Css (0 — '^s expect the large values of 

Yj: to be those that involve the near-zero modes with 
their large-chirality values. Indeed 



E 

T 



(53) 



This expectation is borne out by the numerical results. 

5 5 

On averaging over |(5| = 1 gauge fields, Y^-^ , F22 ' ^'^d 

5 

Y^2 s-re all equal, being the "typical" product of over- 
laps for two near-zero modes. (Yj^ )|Q|=i 1^ shown as a 
function of T in Fig. [9| Results for modes —1 and —2 
from the near-zero-mode quartet match these because, as 
discussed above, the chirality of mode —1 is identical to 
that of 1 and —2 to that of 2. Thus, the sum over all 

the zero modes, i,j e {±1, ±2}, of Y^j gives 4 x 4 = 16 
times the square of the chirality for a typical zero mode. 
This is divided by 16 in the contribution to the discon- 
nected correlator, as in Eq. (49), and so the contribution 



becomes exactly what is required to match that from the 
one zero mode for continuum quarks, up to a renormal- 
ization factor for the taste-singlet pseudoscalar current. 

The nonzero modes, for example mode 3, have small 
chirality and therefore Y^^ is small, as also shown in 
Fig. [9] In the continuum this would be zero. Here it is 
not zero for nonzero lattice spacing but tends to zero as 
a — )• 0. In fact, because we find that Yrr — Kjs ~ 
for modes in the nonzero-mode quartet, r, s € {3, 4, 5, 6}, 
then the total contribution from the quartet, when di- 
vided by 16, cancels against the contribution from X^^. 
divided by 4 in the total pseudoscalar flavor-singlet cor- 
relator, as for the near-zero-mode quartet. 

It is also worth discussing the cross-term Y^^ between 
the near-zero-mode quartet and the nonzero-mode quar- 
tet since this would also be identically zero in the con- 
tinuum. Figure [9] shows the results for Y^^ (T), which, 

summed over T, has a value which is the square root 

5 5 

of the product of the sums over T of F^i and Y^^ . The 
lower plot of Fig.|9]then shows explicitly how J2t ^13 (^) 
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vanishes as a ^ 0. Similar behavior is seen for other 
terms that are related to the chirality of nonzero modes. 

For X'^^ the results for the diagonal case are the same 
as for y . The results for the off-diagonal X"' are less 
clear a priori. In fact, we find in all cases that the off- 
diagonal correlated overlaps within a quartet are zero 



when averaged over gauge fields. Figure 10 illustrates this 

for modes 1 and 2 in the near-zero-mode quartet. X^-^ 

and X22 are large (being equal to Y^^ and F22 ) but XJ2 
has an average of zero. The same results are obtained for 
the —1 and —2 modes. We also see an average of zero for 
the correlated overlaps between the positive and negative 
eigenmodes within the quartet. This is illustrated for 



-1 in Fig 
The size of 



10 



X-^ _2 is very similar. 



' correlated overlaps changes very little with 
the lattice spacing. Figure 10 also compares correlated 



overlaps X^^ on the coarsest lattices, set 1, with those 
for the finest lattices, set 5. 

An average of zero is also seen for off-diagonal terms 
between modes in the first nonzero-mode quartet and be- 
tween modes in the near-zero-mode quartet and modes 
in the first nonzero- mode quartet. These points are illus- 



trated in Fig. 11 



To understand X^^ more completely, we must also 
study correlated overlaps between positive and negative 
nonzero eigenmodes. Although this is not relevant to the 
behavior of the 't Hooft vertex, it shows very clearly how 
the connected contribution to the r/ correlator becomes 
equal to that of the tt meson in the continuum limit, up 
to a renormalization factor that arises because the taste- 
singlet pseudoscalar current is not absolutely normalized. 
Figure 



12 



shows the correlated overlaps, X^^ , between 
the mode r = 3 and all the negative modes that cor- 
respond to the first negative nonzero quartet (which is 
the Ue "mirror" of the first positive nonzero quartet), 

i.e., s = —3, —4, —5, and —6. Interestingly, the corre- 

5 5 

latcd overlaps that arc nonzero here are ^3 ^5 and X^ _q. 
These are equal and each about half the size of X^ ^ = 

5 , _ r~~i 

and Fig 



10 i 
tor r 



Y^i (compare Fig 
correlated overlap 
having the same size as X^ 



12). Likewise, the nonzero 
5 appears with s — —3 and —4 
-5, —6. The correlated 



overlaps Xjg and Xq^ show the corresponding pattern. 
Note the parallel with what happens in the case of the 
near-zero-mode quartet, where ^7-i ^^d X^ 
zero (as shown in Fig. 



10) 



and ATg _ 



2 tend to 



3 and ^3 _4 tend to 



zero too. The differencenere is that another pair belongs 
to the mirror quartet, whereas 1, 2, —1, and —2 form a 
single quartet that is its own mirror. The pattern seen 
in Fig. [12] is repeated for other nonzero-mode quartets. 



0.014, while Xl_j and 



For example, X^ _g and X^_ 

^7 -s much much smaller. 

Thus, the large contributions from nonzero modes 
to the connected correlator for the taste-singlet pseu- 
doscalar meson come from correlated overlaps connecting 
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FIG. 13: A histogram of values of the combination 
— X(Tniid)/4 + y(Tjnid)/16 calculated from the near-zero 
modes i,j = ±1,±2 for the \Q\ — 1 configurations for sets 
1, 3, and 5. The results are plotted for time separation, Tmid, 
set to the midpoint of the lattice. 



members of a quartet and members of its mirror quartet, 
in fact members of the opposite pair of the mirror quar- 
tet. When these correlated overlaps are summed over a 
quartet they give a result, per quartet member, approxi- 
mately equal to that of a typical near-zero-mode contri- 
bution. The near-zero-mode contributions, on the other 
hand, come from diagonal terms, as a result of nonzero 
chirality. On adding all modes together, as in Eq. (47 1, 



and dividing by 4 we obtain a result per quartet, similar 
to that in Eqs. (50) and (52). The way in which this is 



achieved is rather different from that for the Goldstone tt 
meson, and the different mode contributions follow more 
closely that of the continuum. A difference with both the 
continuum and the staggered Goldstone tt is that there is 
a constant of proportionality which is the square of the 
chirality of the zero modes. The disconnected terms can- 
cel all diagonal connected contributions (having in fact 
the same constant of proportionality) , and therefore we 
finally obtain, for the 77' correlator, a result that tends 
to Eq. (51 ) in the continuum limit, once the taste-singlet 



pseudoscalar current is appropriately normalized. 

Let us now demonstrate the cancellation between the 
connected and disconnected contributions from the near- 
zero modes more explicitly. Figure [13] shows histograms 
in the |(5| = 1 sector for 



E 

ij = ±l,±2 



X^, (T) V, (T) 
4 16 ^ 



(54) 



evaluated at the midpoint of the lattice, T = Tmid for 
the three sets — 1, 3, and 5 — that have the same physical 
volume but different lattice spacings. Then Tmid corre- 
sponds approximately to the same physical time sepa- 
ration in each case. From Fig. [13] it is clear that this 
combination of X and Y, which skeptics have worried 
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FIG. 14: Overlaps X^J averaged over the 27 lattices from the 
finest ensemble, set 5, that had topological charge \Q\ — 2. 
Results are shown for mode 1 from the first near-zero-mode 
quartet and mode 3 from the second near-zero-mode quartet. 



could be troublesome, is in fact zero on average at ev- 
ery value of the lattice spacing. The histogram of values 
shows that the distribution is somewhat broader on the 
coarser lattices, but there is no other effect from the lat- 
tice spacing. 

In the above discussion, we have focused on the |(5| = 1 
case because that is the easiest one with which to study 
near-zero and nonzero modes. However, results on con- 
figurations with other Q values also behave exactly as ex- 



pected from this picture. Figure 14 shows results for cor- 



related overlaps X^j for 27 configurations with |(5| = 2 
from the finest, set 5 lattices. The correlated overlaps are 
between modes 1 and 3 which are now members of two 
separate near-zero-mode quartets. We see that there is 
negligible correlated overlap between modes from differ- 
ent near-zero quartets, so the counting for each quartet, 
taken care of by the subsequent division by 4 for the con- 
nected contribution, is exactly as for the IQI = 1 case. 



2. Flavor- singlet scalar mesons 

The flavor-singlet scalar case is easy to analyze both 
in the continuum and for staggered fermions because of 
the simple form of the taste-singlet scalar, F = 1/. The 
orthogonality and normalization of the eigenvectors give 

X]tCrs(0 ^ ^rs- Thus, the disconnected contribution in 
the continuum becomes: 



{i\r + m,){iXs + m) 



(55) 



where the sum is over all eigenmodes. The connected 
contribution is 



1 



[iXr + m)2 



(56) 



FIG. 15: Diagonal and off-diagonal correlated overlaps (T) 
between eigenvectors 1 and 3 that contribute to the discon- 
nected piece of the flavor-singlet scalar meson correlator at 
zero spatial momentum. Results are given for the average 
over IQI = 1 configurations in the fine ensemble, set 5. 



and we see that it is canceled by diagonal terms from 
Eq. (55). In particular, for |Q| — 1, the single zero- 
mode contribution to the total flavor-singlet correlator 
cancels between V and C to give a finite result for A^*^ as 
TO 0. For staggered fermions, Eqs. (55) and (56) still 



hold, with a sum over the total number of eigenmodes. 
By taking a suitable average over the eigenvalues in a 
quartet, the cancellation of diagonal terms quartet by 
quartet mimics that of the continuum. In particular, 
neglecting the near-zero relative to to and dividing 
C by 4 and V by 16, it is clear that exactly the same 
cancellation of the contributions from the near-zero-mode 
quartet occurs as in the continuum. 
Figures 



15 



and 



16 show a representative sample of YK 



15 



and Xlg, plotted as a function of T for set 5. Figure 
shows correlated overlaps Yj^ for a near-zero mode and 
Ygg for a nonzero mode as well as the off-diagonal F/g. 
Set 5 lattices have a time extent of 20, so we expect values 
around 0.05, such that the sum over T yields 1. The re- 
sults for all members of the near-zero-mode quartet agree 
with those of Yj\ and those of the first nonzero-mode 
quartet agree with those of ¥33 . Unlike the pseudoscalar 
case, is not zero but as large as Y]\ and 1^33 since 
X^T ^rsi"^) ~ 1 ^^or all r, s both in the continuum and on 
the lattice. 



Figure 16 shows the correlated overlaps Xlg{T). Here 
-^rsi^)— ^rs, and it is clear that the diagonal cor- 
related overlaps are the same as those of the appropriate 
Y^^ and the off-diagonal correlated overlaps are zero in 
each case. 

which gives the Xl^ 
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Further detail is shown in Fig. 
between nonzero modes in mirror quartets, and between 
the positive and negative eigenmodes of the zero mode 
quartet. Some of these correlated overlaps are large in 
the pseudoscalar case. None of them is large here and all 
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FIG. 16: The correlated overlaps X^g{T) between near-zero 
modes 1, 2 and —1 (top) and between near-zero mode 1 and 
nonzero mode 3 (bottom) that contribute to the connected 
piece of the flavor-singlet scalar meson correlator at zero spa- 
tial momentum. Results are given for the average over |Q| = 1 
configurations in the fine ensemble, set 5. 
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FIG. 17: The correlated overlaps X^^{T) between nonzero 
mode 3 and modes —3, —4, —5, and —6 from its mirror quar- 
tet. Results are averaged over \Q\ = 1 configurations for set 5. 
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FIG. 18: A histogram of values of the combination 
-X^(rmid)/4 + y^(Tmid)/16 calculated from the near-zero 
modes for the |Q| = 1 configurations for sets 1, 3, and 5. The 
results are plotted for time separation set to the midpoint of 
the lattice, Tmid- 



yield zero after summing over T. Quite different behavior 
is seen in the different correlated overlaps, however. In 
particular, we see once again in these correlated overlaps 
the distinction between different pairs in the nonzero- 
mode quartets. 

Figure 18 shows histograms in the |Q| = 1 sector for 



E 

=±1,±2 



xi.(r) yA(r) 

'4 + 16 ' 



(57) 



evaluated at Tmid for the three sets — 1, 3, and 5 — that 
have the same physical volume but different lattice spac- 
ings. From Fig. [18] it is clear that, as in the pseudoscalar 
case, this combination of X and Y , which corresponds 
to the potentially divergent contribution of the near-zero 
modes to the scalar meson correlator, again actually van- 
ishes on average at every value of the lattice spacing. 
The width of the histogram distribution is the quantity 
which changes with lattice spacing, becoming more nar- 
rowly peaked around zero as the lattice spacing goes to 
zero. 



3. Flavor-smglet vector, axial vector, and tensor mesons 

The correlated overlaps for the flavor-singlet tensor 
case behave similarly to the pseudoscalar and scalar. No 
simple analysis of correlated overlaps in terms of the 
chirality or normalization of the modes is possible and, 
indeed, we find that none of the correlated overlaps is 
large. Figure [19] shows that the key requirement for a 
sensible flavor-singlet correlator holds, i.e., that the off- 
diagonal correlated overlaps between different members 
of the near-zero-mode quartet are consistent with zero. 
This means, as above, that the connected and discon- 
nected near-zero-mode contributions cancel rather than 
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FIG. 19: The correlated overlaps X^j (T) between near-zero 
modes 1 and 2 that contribute to the connected piece of the 
flavor-singlet tensor meson correlator at zero spatial momen- 
tum. Results are given for the average over = 1 configu- 
rations in the fine ensemble, set 5. 



FIG. 20: The correlated overlaps FJ (T) for fi — x between 
near- zero modes 1 and —1 that contribute to the disconnected 
piece of the fiavor-singlet vector meson correlator at zero spa- 
tial momentum. Results are given for the average over \Q\ = 1 
configurations in the fine ensemble, set 5. 



x'^, ■ — X — i x'^ 



giving a potentially divergent piece. 

The flavor-singlet vector and axial-vector cases behave 
somewhat differently, which can be traced back to the 
fact that the taste-singlet versions of these operators cou- 
ple even and odd lattice sites together rather than even- 
to-even or odd-to-odd as with the other examples. The 
axial vector and vector behave in the same way, so we 
only show results here for the vector case. As discussed 



0.0001 



in Sec. Ill in the continuum there is no zero-mode con- 



tribution to the disconnected piece of the flavor-singlet 
vector meson correlator because and 7^ anticommute. 
We show in Fig.[20|how this works for staggered fermions. 

-7" 
Because d 



couples odd and even sites, and /_i has the 

opposite sign on odd sites to /i, then i^^-^ and C-1-1 have 
opposite sign. This means that the near-zero-mode con- 
tribution to the disconnected correlator from Y7i has 
opposite sign to that from yjLi- This is seen clearly 
for fj, = X in Fig. 20 Summing over z,j G {±1,±2} 



then clearly gives a total disconnected contribution to 
the flavor-singlet vector meson correlator of zero. 

The connected contributions are not zero time slice 
by time slice, as we see from Fig. [21] Correlated over- 
laps X^^ (equal to their counterparts in Fig. 20) 
become pure oscillations, (—1)'^, that cancel upon sum- 
ming over T. Oscillatory terms are a feature of staggered 
meson correlators, stemming from opposite-parity contri- 
butions to the correlator and, ultimately, the remaining 
time doubling. They do not then affect the properties 
of the ground state meson, in this case the flavor-singlet 
vector meson. The off-diagonal X^j are close to zero 
and also oscillatory. Thus, once again there is no signifi- 
cant net contribution from near-zero modes to the flavor- 
singlet vector meson. 




FIG. 21: The correlated overlaps X^J (T) for jj, = x between 
near-zero modes 1 and 2 that contribute to the disconnected 
piece of the fiavor-singlet vector meson correlator at zero spa- 
tial momentum. Results are given for the average over = 1 
configurations in the fine ensemble, set 5. 



4. Summary 

We conclude that the behavior of the staggered flavor- 
taste singlet meson correlators in every case follows that 
expected in continuum QCD. In particular, no chiral- 
limit divergence results from the near-zero modes. In all 
cases, we find that the take the same average value 
for i,j both in a given quartet. Thus, the disconnected 
contribution from a quartet of degenerate eigenvalues is 
4 X 4 = 16 times that of a single mode. The off-diagonal 
correlated overlaps Xij are zero for i,j within a quartet in 
every case. This means that the connected contributions 
give instead 4 times that of a single mode. Including the 
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factors of 4 and 16 in Eq. (49) means that the correla- 



tor is effectively made of single-species contributions, as 
in continuum QCD, and, in particular, the contribution 
from near-zero modes cancels as it does there. 

The pseudoscalar and scalar cases are particularly sim- 
ple to analyze, both for zero and nonzero modes, and to 
see the clear correspondence with continuum behavior. 
The correlated overlaps for the taste-singlet pseudoscalar 
between mirror quartets are a striking demonstration of 
how staggered fermions conspire to give the "right" an- 
swer, but sometimes in a rather nontrivial way. The 
match demonstrated between the taste-singlet and the 
Goldstone pseudoscalar also leads to a practical sugges- 
tion that may improve the determination of the rj mass 
using staggered fermions. The calculation of the taste- 
singlet connected and disconnected contributions is par- 
ticularly statistically noisy because of the point-split na- 
ture of the taste-singlet operator. It may be preferable, 
although numerically challenging, to determine instead 
the near-zero-mode eigenvectors and then subtract their 
contribution from the Goldstone pseudoscalar correlator. 
This must agree with the 77' correlator in the continuum 
limit and yet is constructed of local operators and so has 
significantly less gauge noise. 



V. CONCLUSIONS 

This paper adds weight to the growing evidence that 
shows that staggered fermions behave in the correct way 
to reproduce QCD in the continuum limit, even with the 
rooted determinant. Here we have focused on the eigen- 
vectors of the staggered-fermion Dirac operator and the 
way in which the 't Hooft vertex and flavor-singlet me- 
son correlators are built from the overlaps between dif- 
ferent eigenvectors, using the appropriate taste singlets. 
The important overlaps are those between eigenvectors 



within a near-zero quartet since these could have gen- 
erated dangerous singular terms as m — >■ 0. From our 
theoretical results we determine a condition for the local 
overlaps that needs to hold and then test this numeri- 
cally and demonstrate that it does. Indeed we see that 
the near-zero-mode quartet in all cases behaves function- 
ally in such a way to reproduce the required behavior of 
four copies of a single mode that mimics the expected 
behavior in the continuum. 

Most of our results are not surprising, but in provid- 
ing a clear link between the theoretical requirements and 
the numerical results for the eigenvector overlaps, we add 
further confidence to the soundness of the framework for 
the accurate phenomenology that is being done with stag- 
gered fermions. We demonstrate most directly that a cal- 
culation of flavor-singlet meson masses, notably that of 
the rj meson, should give the correct QCD result. This 
is not a substitute for doing the full calculation and this 
is underway [SI] |32] . 
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Appendix A: Detailed Formulas 



For Sec. IV C it is convenient to spread staggered-fermion bilinears over a hypercube, inserting the lattice gauge 
field to preserve gauge invariance. An explicit construction is 



^^(^) = ^IIx(a; + 6)x(a; + 6), 

h 

16 ^-^ 

b 



A'ji^) = -^y^ 'n^'{x + d-b^^''^)-qi{x + b)ri2{x + b)ri'i{x + b)rii{x + b) x 
16 ^-^ 



X{x + d~ U^'^)U{x + d- U^'\x + b)x{x + b), 
Pi{x) = h.^m{x + b)rj2{x + b)T]3{x + b)i]4,{x + b) x{x + d - b)U{x + d-b,x + b)x{x + 6), 



(Al) 
(A2) 
(A3) 



(A4) 
(A5) 
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where b 



but b 



runs over the 
'p — 



-site hypercube with origin x; 



bp, p ^ \x,v\ and d = (1 



= a-^A, A = /I, t/, 
2 + 3 + 4) a. Gauge 



invariance is ensured via averages of parallel transport 
over paths from x to x' , U{x, x').^ Under shift symmetry 
these are all taste singlets. The vector current and scalar 
density satisfy the Ward identity corresponding to quark- 
number conservation for all a, and the axial-vector and 
pseudoscalar density satisfy the anomalous Ward identity 
as a — >■ 0. In practice, we use in place of U the HISQ- 



smeared gauge field W, defined in Eq. (B5) below 



For brevity and clarity, it is then helpful to write 



Si{x) 




(A6) 


vrix) 




(A7) 




- - LLU 

= X^(^I X, 


(A8) 


A^iix) 


- - u5 

= *X7> X, 


(A9) 


Pi{x) 


= X7/X, 


(AlO) 



which with Eqs. (|Al|)-(|A5|) define 1/ 

v5 



7/ 



ll 



and 



7j, when acting on x, Xj an d the eigenvectors of Dgtag 



Sec. 



IV C 



for the analysis of m 

When constructing the correlators for the rf and other 
flavor-taste-singict correlators, it is more customary to 
restrict the operators to one time slice. In Sec. |IV D 1[ 
therefore, we average over spatial cubes only [and then 
time slices, cf. Eqs. (46)]. Vi and Aj remain as in 
Eqs. (A2) and (A4|, because they naturally extend over 
a timelike link or three-dimensional cube. For P/, the 
operator is defined as attached to point on a time slice 
by averaging over all hypercubes that have a corner at 
that point, i.e., extending both forwards and backwards 
in time. Then averaging over a time slice is straightfor- 
ward. 



Appendix B: Improved Staggered Actions 

To introduce improved staggered-fermion actions, it 
is convenient to proceed in steps, introducing notation 
along the way. The first step is "Fat 7" smearing ^ , 



sym 



(Bl) 



which yields paths of length 3, 5, and 7. Here 
T±pUi_,{x) = U±p{x)Up,{x ± pa)U^p{x ± pa), U-p{x) = 
Uj,{x-pa). It is easy to check that the smearing intro- 
duces a form factor that reduces the coupling to taste- 
changing gluons [32] . 

As is often the case with smearing algorithms, Fat7 
smearing introduces additional discretization errors. 



^ There is no bar on U{x + fe^^'^x + b) because only the one-link 
path enters. 



These can be removed by introducing an order-a^ im- 
provement 73' 



where 



Vp = {Tp - \Cp) Up 



CpUp^Y.^Tp-T^pfUp, 



(B2) 



(B3) 



introduces the five-link Lepage term. The discretization 
error of the simple difference operator in Eq. ^ can be 
removed with the three-link Naik term |70j . 

5'Naik = -^a^^-qp{x)x{x) (Tp - T^pf x{x), (B4) 

where now T±px{x) — U±p{x)x{x ± /ia). 

For the HISQ action, Fat7 smearing is applied twice, 
with the Lepage correction taken at the second step 



Wp = {Fp-\Cp)UFpUp, 



(B5) 



where U denotes a reunitarization and projection to 
SU(3). [The SU(3) projection makes little difference in 
practice.] The HISQ action is then 



'S'hISQ — Sstag^{Wp) + S^s,i]^(plFpUp), 



(B6) 



substituting for the original gauge field Up as shown. 

For completeness we write the Fat7xAsqtad [30] and 
Asqtad [73] actions in this notation: 



5, 



Fat7xA sqtad S-stag {Wp)+ ^Naik {UTp Up ) , (B7) 

Wp^ {Fp-\Lp)UFpUp, (B8) 

Sa sqtad — 5'stag(V"^i) + S^a.ik{Up), (B9) 



Unfortunately, Ref. [27] referred to Fat7x Asqtad as 
"HISQ." The Asqtad action defines the rooted determi- 
nant in the MILC ensembles [TT] [5S], which have been 
used by the zero-temperature results cited in the Intro- 
duction. For this action there is an additional tadpole- 
improvement step in which one replaces T±pX and T±pUp 
by Uq^T±pX and UQ^T±pUp, respectively, where uq is a 
measure of the mean link. In MILC's simulations of the 
Asqtad action [TT], ug is set by the fourth root of the 
1x1 Wilson loop (the plaquette). (The reunitarization 
in HISQ makes tadpole improvement unnecessary.) 



Appendix C: Further remarks on Refs. [35', '36] 

Creutz |351 136] makes several remarks that sound sim- 
ple, and thus seem to be accepted by nonexperts, but 
they do not withstand careful scrutiny. One, explained 
elsewhere [TO], is that the different tastes have differ- 
ent chirality. As discussed above, all near-zero modes 
within a common quartet possess (identically for mir- 
rors; empirically otherwise) the same taste-singlet chi- 
rality, Eqs. (25 1 or (26). The nonzero modes all have 
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(nearly) zero taste-singlet chirality. Finally, all modes 
have no net Goldstone chirality, J2x fsi^)lpfsix) — 
e{x)fs{x)Mx) = f -six) f six) = 0. 
Another incorrect statement [36 concerns the an- 
gle of the strong CP problem, which can appear via a 
modified mass term 

mijjip 1-^ mcos{9) ipip + imsm{9) ^pj^ip. (CI) 

Creutz states correctly that 9 obtains a physical meaning 
via the anomaly and, hence, the ultraviolet regulator. He 
also states, incorrectly, that staggered fermions cannot 
possess this property, owing to the exact symmetry. 

With this symmetry, the following two mass terms are, 
of course, equivalent: 

mxxix) ^ mcos{ip) xxix) + 'i'msmi(p) eix)xxix)- (C2) 

In the continuum limit, however, this corresponds to 

mqqix) O m cosi(p) qqix) + im sin{ip) qj^^^^qix) , (C3) 



namely is a taste nonsinglet. It is superficially the 
kind of transformation used to set up twisted-mass 
Wilson fermions |76| . but we have seen no argument 
that proves it is the same. In particular, unlike the 
unsubstantiated mass terms posited in Refs. ^521 1^ . 

^i^) — fsix)f-si^) is off-diagonal in any basis where 
Eq. (fTsFmakes sense. 



The correct analog of Eq. (CI I is 



mxxix) ^ mx [cos(6') -f i sin(6')7|] xix)- (C4) 



The taste singlet 7I extends across a hypercube and de- 
pends on the lattice gauge field. It thus relies on the 
regulator for its definition, as it must. To simulate the 
9 vacuum via the fermion mass, one needs to implement 



Eq. (C4|, not Eq. (C2) [251E2] 
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